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ABSTRACT 

Observations of damped Lya absorbers (DLAs) can be used to measure gas-phase metallicities at 
large cosmological lookback times with high precision. Furthermore, relative abundances can still 
be measured accurately deep into the reionization epoch (z > 6) using transitions redward of Lya, 
even though Gunn-Peterson absorption precludes measurement of neutral hydrogen. In this paper 
we study the chemical evolution of DLAs using a model for the coupled evolution of galaxies and 
the intergalactic medium (IGM), which is constrained by a variety of observations. Our goal is to 
explore the influence of Population HI stars on the abundance patterns of DLAs to determine the 
degree to which abundance measurements can discriminate between different Population III stellar 
initial mass functions (IMFs). We include effects such as inflows onto galaxies due to cosmological 
accretion and outflows from galaxies due to supernova feedback. A distinct feature of our model is 
that it self-consistently calculates the effect of Population HI star formation on the reionization of an 
inhomogeneous IGM, thus allowing us to calculate the thermal evolution of the IGM and implement 
photoionization feedback on low-mass galaxy formation. We find that if the critical metallicity of 
Population HI to II/I transition is < 10~*Zq, then the cosmic Population HI star formation rate (SFR) 
drops to zero for z < 8. Nevertheless, at high redshift (z ^ 6) chemical signatures of Population HI 
stars are retained in low mass galaxies (halo mass < 10^ M0). This is because photoionization 
feedback suppresses star formation in these galaxies until relatively low redshift (z ~ 10), and the 
chemical record of their initial generation of Population HI stars is retained. We model DLAs as 
these low-mass galaxies, and assign them a mass-dependent H i absorption cross section in order 
to predict the expected distribution of DLA abundance ratios. We find that these distributions are 
anchored towards abundance ratios set by Population II supernovae yields, but they exhibit a tail 
which depends significantly on the Population HI IMF for z > 5. Thus, a sample of DLA metallicity 
and relative abundance measurements at liigh-rcdshift holds the promise to constrain Population HI 
enrichment and the Population HI IMF. We find that sample of just 10 DLAs with relative abundances 
measured to an accuracy of 0.1 dex is sufhcient to constrain the Population III IMF at 4(7. These 
constraints may prove stronger than other probes of Population HI enrichment, such as metal-poor 
stars and individual metal-poor DLAs. Our results provide a global picture of the thermal, ionization, 
and chemical evolution of the Universe, and have the potential to rule out certain Population HI 
scenarios. 

Subject headings: Cosmology: dark ages, reionization, first stars - galaxies: abundances - galaxies: 
evolution - galaxies: ISM - galaxies: quasars: absorption lines - stars: Popula- 
tion HI 



1. INTRODUCTION 

The first stars in the Universe formed out of the pri- 
mordial interstellar media (ISM) of the first galaxies and 
are therefore expected to be metal- free. Formation of 
these so-called Population HI stars is an important mile- 
stone in cosmic evolution as the UV photons produced by 
these objects presumably initiated the process of reion- 
ization of the intergalactic medium (IGM). Reionization, 
and the consequent reheating of the IGM result in new 
physical processes such as photoionization feedback that 
are expected to influence subsequent galaxy formation 
and evolution. Apart from reionization. Population HI 
stars also contribute to chemical feedback. The metal 
yield of these stars pollutes the ISM and creates condi- 
tions for formation of metal-enriched Population II stars 
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that we observe today. Population HI stars could also 
seed the formation of supermassive black holes seen in 
centres of galaxies at low redshift. These stars thus 
form a significant piece of the puzzle of galaxy forma- 
tion (Bromm & Larson 2004; Loeb 2010). 

The initial conditions of galaxy formation arc 
well known thanks to cosmological constraints at 
the last scattering surface (Ciardi & Ferrara 2005; 
Bromm & Yoshida 2011). In the first galaxies that result 
from these initial conditions (at z > 20 in dark matter 
haloes with mass of about 10^ M©), gas can fragment to 
form Population HI stars. The chemistry of gas in these 
haloes is simple, and formation of the first generation of 
stars is unaffected by complicated effects like feedback. 
Physical factors such as magnetic fields are expected to 
play a negligible role in the formation of Population HI 
stars. We can therefore hope that the physics of Pop- 
ulation HI star formation is simpler than that of star 
formation in the present day universe. 

Nonetheless, properties of Population HI stars, such as 
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their mass, are not well-understood, primarily because 
of the wide range of scales involved in the fragmenta- 
tion and collapse of gas in high redshift haloes. On the 
galactic scale, the mass of the first galaxies is debatable 
(Bromm & Yoshida 2011). Since the molecular hydrogen 
cooling threshold provides a lower bound on halo mass for 
gas cooling and fragmentation, we expect that the first 
stars will be hosted by haloes with virial temperature 
Tvir = 10'^ K. However, molecular hydrogen is fragile and 
can be easily destroyed as the IGM is transparent at the 
required energies. If the molecular hydrogen is destroyed 
by the Lyman- Werner background produced by the very 
first few stars, then the formation of new stars could be 
delayed up to a lower redshift in all but the most massive 
haloes (e.g., Barkana & Loeb 2001). On the proto-stellar 
scale, it is unclear whether the gas collected in these 
haloes will fragment, or what will be the mass of the 
fragments if it does. It has been suggested that dynami- 
cal effects can cause fragmentation even when metal-line 
cooling is not available (Dopcke et al. 2012). But the 
magneto-hydrodynamics in this regime is very difficult 
to simulate. In summary, both the initial mass function 
(IMF) as well as the cosmic star formation rate (SFR) 
history of Population III stars remain unclear. We are 
thus forced to ask if there are any empirical constraints 
on these quantities. 

Unfortunately, no understanding of Population III 
stars has emerged from the observational side cither. 
Currently, two lines of observational inquiry have aimed 
at an unambiguous detection of the chemical signatures 
of Population III stars. These either probe (1) low mass 
extremely metal-poor stars in the Galaxy's halo, or (2) 
metal-poor neutral gas reservoirs, presumably galaxies, 
seen in quasar absorption spectra at high redshift (z > 2) 
as damped Lya absorbers (DLAs). In the first case. 
Galactic stellar archaeology has been considered as an 
important probe of high redshift star formation envi- 
ronment (Frebel 2011). These long-lived stars could 
potentially carry signatures of enrichment by Popula- 
tion III supcrnovae (Rollinde et al. 2009). However, no 
such detection has been made, although current obser- 
vations have a potential sensitivity of metallicity down 
to [Fe/H] = -9.8 (Frebel & Norris 2011).^ Moreover, 
although they are a valuable tool for constraining Galac- 
tic chemical evolution, it is not clear if metal-poor stars 
are an unambiguous probe of the cosmic SFR and IMF 
of Population HI stars, because these observations are 
limited to the peculiar assembly history of the Galaxy 
and because abundance ratios in a single metal-poor 
star cannot probe the whole Population HI IMF. In the 
second case, metal-poor DLAs at z ~ 2 likely repre- 
sent systems that form in peculiar, pristine patches of 
the IGM. They could be direct descendants of high red- 
shift haloes and, due to inefficient star formation, prob- 
ably retain signature of Population HI star formation 
(Salvadori & Ferrara2009, 2012). However, observations 
of such systems have also failed to show up chemical 
abundance patterns that cannot be explained with Popu- 
lation II yields alone (Cooke et al. 2011; Fumagalli et al. 
2012). 

In the absence of current observational constraints, 

3 The lowest-metallicity star known has [Fe/H] = -4.88 ± 0.12 
when computed using a 3D model atmosphere(CafTau et al. 2011). 



there are suggestions that future facilities such as the 
James Webb Space Telescope (JWST) and the Large 
Synoptic Survey Telescope (LSST) may detect pair- 
instability supernovae (PISN) from any high-mass Pop- 
ulation HI stars in situ at high redshift (Pan et al. 
2012a,b). It has also been argued that massive Popula- 
tion HI stars could trigger long GRBs (Bromm & Loeb 
2006). Detection of these objects will also potentially 
probe the first stars (Levesque 2012). 

In this paper, we study the possibility of probing the 
Population HI SFR and IMF by measuring chemical 
abundance ratios in low mass galaxies at high redshift 
{z ~ 6). We set up a conservative model of evolution 
of galaxies and the IGM, constrained by low-redshift ob- 
servations, and show that different Population HI IMFs 
leave distinct signatures in the metal abundance ratios of 
low mass galaxies at high redshift. Further, these galax- 
ies can be observed as DLAs in the spectra of any back- 
ground quasars. These absorption-line observations can 
be used to measure gas-phase metallicities in these galax- 
ies with high precision. Even deep into the reionization 
epoch (z > 6), abundance ratios can still be measured 
accurately in such systems using transitions redward of 
Lya. We model DLAs as these low-mass galaxies by 
assigning them a mass-dependent H i absorption cross 
section, and predict the expected distribution of DLA 
abundance ratios at various redshifts. We find that at 
high redshift the distribution of high redshift DLAs in 
the abundance-ratio space can be uniquely mapped to a 
combination of Population III SFR and IMF. This can 
be a powerful probe of the first stars. 

We describe our modelling of the evolution of galaxies 
and the IGM in Section 2. Our results are presented 
in Section 3. We discuss some caveats in Section 4 and 
present a summary of conclusions in Section 5. 

2. MODELLING GALAXY AND IGM EVOLUTION 

Given a dark matter halo of mass Mh at z = 0, we can 
calculate its complete mean assembly history. This can 
be done using the extended Press-Schechter formalism, 
or using large volume N-body simulations of structure 
formation. We prefer the latter approach for its accu- 
racy and use fitting functions provided by Fakhouri et al. 
(2010). This lets us calculate the mean dark matter mass 
M of a halo at any stage of its assembly history. We 
use this to set up a model of baryonic evolution inside 
the halo. In simplest terms, this model assumes that 
a halo (1) accretes baryons through cosmological accre- 
tion, (2) forms stars from any gas contained in the halo 
for sufficiently long duration, and (3) ejects baryons due 
to supernova-induced outflows. 

With these assumptions, we can associate with each 
dark matter halo a stellar mass and a gas mass Mg 
at any point in its assembly history. Further, using stel- 
lar population synthesis, we can also calculate the total 
amount of metals Mz and the rate of UV photon pro- 
duction Tiph in a halo at any stage of its assembly history. 
All of these quantities are functions of the redshift-zero 
halo mass {Mh). We repeat this exercise for a num- 
ber of logarithmically spaced halo masses between Mj^°^ 

and Af^'^*^, two values that are chosen using a conver- 
gence criterion. The integrated UV photon yield of all 
halos in this mass range is then used in a reionization 
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model to calculate the ionization and thermal history of 
the IGM. Similarly, the integrated outflow yield is used 
with a chemical enrichment model to calculate the aver- 
age IGM mctallicity. The IGM ionization, thermal, and 
chemical evolution models feed back into the galaxy for- 
mation model through photoionization feedback, which 
controls the thermal cut-off on the mass of star-forming 
haloes, and by changing the metallicity of inflows from 
the IGM to the galaxy. 

This galaxy formation picture is often referred to as 
the "bathtub" model, and has been used to study various 
problems in galaxy formation to a reasonable degree of 
accuracy (Bouche et al. 2010; Krumholz & Dekel 2012). 
We will now describe our method in detail. A summary 
of free parameters of our model is presented in Table 2. 

2.1. Halo growth 

The mean mass growth rate of dark matter haloes in 
the ACDM cosmology can be parameterised by the fitting 
function (Fakhouri et al. 2010) 



A/(z) = 46.1^(1 -fl.lz) 
yr 

X ^f7™(H-z)3 + 17A 



M{z) 



(1) 



where M(z) denotes the mean mass in the main branch 
of the merger tree of a halo at redshift z. This growth 
rate corresponds to the mean rate of growth of halo 
mass. This relation can be understood in terms of the 
extended Press-Schcchtcr picture of structure formation 
(Neistein & Dekel 2008). By comparing it with cosmo- 
logical simulations, Fakhouri et al. (2010) have shown 
that it is accurate over a mass range spanning at least five 
orders of mag nitude (10^° M© to 10^^ M© at z = 0) and 
a wide range of redshifts (0 < z < 15). Integrating Equa- 
tion (1) with the boundary condition M{z = 0) = Alh 
gives the mean mass of a halo at any redshift. Note that 
this describes the evolution of the dark matter mass of 
the halo. 

2.2. Baryons 

We assume, as suggested by hydrodynamical simula- 
tions, that infiows, outflows, and star formation in galax- 
ies are in equilibrium (Dave et al. 2012). We now de- 
scribe these processes and calculate the mean evolution 
of gas mass {Mg), stellar mass (M*), and metal mass 
(Mz) in each halo mass bin. The metal mass term de- 
scribes the mass of the metals in the gas phase (ISM); it 
excludes the metals that are locked up inside stars. 

2.2.1. Gas 

Gas infiows replenish the gas content of haloes as they 
evolve. We assume that the mean gas infiow rate is pro- 
portional to the mean halo growth rate 



A/,,in(^) - /s,in ( ^ ) M{Z 



(2) 



We follow Bouche et al. (2010) and take the constant 
of proportionality to be /g.in = 0.7 , except when the 
halo mass is below the filtering mass Mj^in (described 
below), in which case there is no gas infiow and we set 



fg,in ~ 0. The choice of these values is motivated by the 
comparison of measured SFR in massive high-redshift 
galaxies with predicted gas accretion rates (Genel et al. 
2008; Dekel et al. 2009; Bauermeister et al. 2010). 

Once inside the halo, some of the accreted gas is lost 
to star formation. The amount of gas lost is 



Mg,,i{z) = -V(Z) 



dm Sim 



) ■ '4>[t{z) — '''{ra)] ■ [m ~ mr{m)], (3) 



where the first term on the right hand side is the star for- 
mation rate (which we describe below), and the second 
term accounts for mass loss from evolving stars via stel- 
lar winds and supernova explosions. In the second term, 
4){rn) is the stellar IMF, mr{m) is the remnant mass left 
by a star with initial mass m, T(rn) is its lifetime, and t{z) 
is the cosmic time. The lower limit m; of the integral is 
such that T{mi) = t{z). The instantaneous recycling ap- 
proximation (IRA) is often used in calculating the mass 
loss from stars (e.g., Krumholz & Dekel 2012). In this 
approximation, stars with initial mass above a certain 
value mg are assumed to die instantaneously (r = 0) 
while those with initial mass less than toq are assumed 
to live forever (r = oo). However, although this is a 
reasonable assumption to make while studying, say, the 
solar neighbourhood, it is not a good approximation at 
high redshift, when stellar lifetimes are comparable to the 
Hubble time. Therefore, we do not use IRA in Equation 
(3). 

Gas is also lost due to outflows resulting out of 
supernova-driven winds. This effect is expected to 
be proportional to the star formation rate and in- 
versely proportional to the depth of the halo potential 
well (Daigne et al. 2006; Erb 2008; Bouche et al. 2010; 
Krumholz & Dekel 2012). We write the mass loss due to 
outflows as 

2e /■"" 

dm {4>{m) 

CSC 'Jrrii 

X ^j[t{z) - r(m)] ■ Ei,in{m)} , (4) 



Mg r,ut{z) — 

V 



where v"^^^ = 2GM/Rvii is the escape velocity of the halo, 
and Ekin{m) is the kinetic energy released by a star of 
mass m. The parameter e is flxed by calculating the 
total baryon fraction in structures, /bstructj at z = 
(Fukugita & Peebles 2004). 

Gombining the contribution of above three processes 
of star formation, infiows, and outflows, the evolution of 
the gas mass of a halo can be written as 

Mgiz) = Mg,in{z) + + Mg,ontiz), (5) 

where each term is redshift- and halo-mass-dependent. 
2.2.2. Stars 

We assume that star formation in a halo tracks the 
total amount of gas inside that halo (Bouche et al. 2010). 
Thus 



Mr, 



dyn 



(6) 



where is the star formation efficiency, and tdyn is the 
galaxy dynamical time. This is a significant simplifi- 
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cation of what is a highly comphcated process. How- 
ever, Equation (6) follows from the empirical Kcnnicutt- 
Schmidt (KS) relation (Krumholz & Thompson 2007) in 
the limit of marginally unstable disks. The KS relation 
does not exhibit any evolution till z ^ 2 (Daddi et al. 
2010). Furthermore, although it is possible that changes 
in Equation (6) will result in modified evolution of star 
formation in some mass bins, we do not expect it to have 
a qualitative effect on our results. 
The galaxy dynamical time can be written as 



t.,. = 2 X lOV 



200 kms- 



(7) 



where R1/2 is galaxy disk scale- length and Vc is the halo 
circular velocity (Bouche et al. 2010). We take the disk 
scale-length to be proportional to the halo virial radius 



/2 



0.05 



A 

OJ 



(8) 



where the halo spin parameter A is taken to be 0.07 
(Krumholz & Dekel 2012). The halo circular velocity is 
given by 



Vr = 23.4 



M 



1/3 



l + Z 



1/6 



10 



1/2 



where 



km/s, (9) 



(10) 



and Ac is the halo overdensity relative to the cosmolog- 
ical critical density, given by 



Ac = IStt^ + 82d - 39(f 



(11) 



for the ACDM cosmology, where d = Of„ — 1 
(Barkana & Loeb 2001). The virial radius is given by 



Rv 



0.784 



M 



108/i-iM, 



•'-'0. 



^1/3 






[f}f„187r2j 



-1/3 



10 



h-^kpc. (12) 



The total stellar mass in the halo evolves as 



dm (j){m) ■ ip[t{z) — T(m)] • [m — 'mr{m)], (13) 



where, as in Equation (3), the second term on the right 
hand side accounts for mass loss from stars. 

2.2.3. Metals 

The initial mass of metals in a halo is zero. Metals are 
produced in stars in the halo and are mixed in the halo 
gas when stars explode as supernovae, and due to stellar 
winds. The total mass of metals AIz in the halo is also 



Table 1 

Stellar IMFs used in this paper 



Model 


Pop. Ill IMF 




Pop. II IMF 


Model 1 
Model 2 
Model 3 


1-100 Mq Salpeter 
35-100 Mq Salpeter 
100-300 Mq Salpeter 


poo 


-100 Mq Salpeter 
-100 Mq Salpeter 
-100 Mq Salpeter 



affected by inflows from the IGM and outflows into it. 
The inflow term is given by 



Mz:miz) = ZiGM{z)Mg^in, 



(14) 



where Zigm is the metal abundance in the IGM, which 
we describe below. 
Metal outflows can be described in similar fashion, as 



Mz.OUtiz) = -Zg{z)Mg^oUt, 



(15) 



where Zg = Mz/Mg is the gas metal abundance in the 
halo. 

Lastly, the effect of star formation can be expressed as 
(Krumholz & Dekel 2012) 



dm(j){m) ■ Tpitiz) — T(m)] 



xmpz{m){l-Q. (16) 



Here, pz(jn) is the mass fraction of a star of initial mass 
m that is converted to metals and ejected. The factor 
1 — C takes into account the fact that not all of the newly 
enriched material is going to mix in the ISM. Some of 
it will be directly ejected out of the halo due to super- 
nova explosions. We follow Krumholz & Dekel (2012) 
and write 

C = exp(-A/,,,i2/A/,et), (17) 

where Mji^i2 is the halo mass in units of 10^^ Mq and we 
set the parameters Miot = 0.3 and Q = 0.9. 

The evolution of the metal mass of a halo can now be 
written as 

Mziz) = MzAz) + Mz,ont{z) + Mzm{z). (18) 
2.3. IMF 

We always take the stellar IMF to have the Salpeter 
form, which is given by 

0(m) = 0om-2-3, (19) 

where the constant 0o is chosen such that 



dmm(j){m) = IMq. 



(20) 



We take mo ~ 0.1 Mq and mi = 100 Mq for Popula- 
tion I/II stars. Population HI stars have different IMF 
in each of our models, as we want to discuss the effect 
of Population HI IMF on various quantities. These are 
shown in Table 1. The stellar IMF in model 3 covers the 
pair instability supernova range, whereas models 1 and 2 
explore the effect of AGB stars and core collapse super- 
novae respectively. Each model is calibrated separately 
to reproduce the cosmic star formation rate (SFR) his- 
tory, the IGM Thomson scattering optical depth to the 
last scattering surface (re = 0.089) (Larson et al. 2010), 
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and the hydrogen photoionization rate in the IGM as 
measured from Lya forest observations. 

Our model requires stellar lifetimes because we do not 
assume instantaneous recycling. We take lifetimes for low 
and intermediate mass stars (0.1 < M < 100 Mq) 
from Maeder & Meynet (1989). Lifetimes for stars with 
higher mass (all of which are Population III) are taken 
from Schaerer (2002). The model also requires chemi- 
cal and UV photon yields of stars with different IMFs. 
Woosley & Weaver (1995) have calculated the chemi- 
cal yields of intermediate mass stars (12 Mq < M < 
40 Mq) with different metallicities {Z = 0, lO"'', 0.01, 
0.1, and 1 times solar). We use their results by inter- 
polating between different metallicities and extrapolat- 
ing beyond the mass range for lower and higher stellar 
masses. Chemical yields of Population III stars (which 
are metal-free and high mass) are taken from the cal- 
culations of Heger & Woosley (2002). Stellar spectra of 
Population I/II stars are calculated using STARBURSt99 
(Leitherer et al. 1999; Vazquez & Leitherer 2005) with 
respective metallicities. Synthetic spectra of Popula- 
tion III stars are taken from Schaerer (2002). 

The transition from Population III star formation to 
Population II star formation in any halo is implemented 
via a critical metallicity, Zcrit-, with our fiducial value as 
Zcrit = 10^4 (Bromm et al. 2001; Bromm & Loeb 2003; 
Frebel et al. 2007). When the ISM metaUicity in a halo 
crosses Zcrit, new stars are formed according to a Pop- 
ulation II IMF. We consider the effect of changing the 
value of Zcrit below. 

2.4. IGM metallicity 

Cosmological inflows bring material from the IGM into 
the ISM of galaxies. The extent to which this dilutes 
the metal content of the ISM depends on the metallic- 
ity of the IGM. We follow the model of Daigne et al. 
(2006) to calculate the average IGM metallicity at any 
redshift. This model divides cosmic baryons into three 
reservoirs: (1) intergalactic medium ("IGM"), (2) inter- 
stellar medium ("ISM"), and (3) stars ("str"). We denote 
the mass densities of these as A/igm, -^^^ism and Mgtr re- 
spectively. Note that Mjsm and Mgtr are mass function 
weighted integrals of Mg and il/* defined for each halo 
mass bin above. 

The IGM mass density evolves according to 



dt 



^ -ab{t)+o{t), 



(21) 



where ab is the average rate of accretion of baryons from 
the IGM on to halos, and o{t) is the average rate of 
outflow of baryons from halos into the IGM. The stellar 
mass density evolves as 



dMst 
dt 



*W-e(t), 



(22) 



where is the cosmic average star formation rate, and 
e is rate of ejection of material of stars into the ISM via 
winds and supernovae, which depends on the stellar IMF. 
Lastly, the ISM mass density evolves as 



dAh 



ISM 



dMi 



IGM 



dt dt dt 

To calculate the IGM metallicity, in addition to the 
total masses, the mass fraction of individual elements in 



each reservoir is also evolved. For an element i, the mass 
fraction in the IGM is Xj^^ = MI™/Migm and that in 
the ISM is Xj^^ = Mj^^/MisM- It can then be shown 
that mass fractions in the IGM evolve as 

jyIGM 

- O (Aj - Aj j /AilGM, 

and those in the ISM evolve as 
dX}^^ 



(24) 



dt 



= [at (Xf ^ - ^) 



+ (e,-eXf^)]/M, 



ISM- 



(25) 

Here, o is the outflow rate from the ISM and is the 
rate at which metal i is ejected by stars. In order to 
completely specify the model, all that is left to specify is 
(1) initial conditions for above equations, and (2) method 
of calculation for the three quantities o, ab, and e. 
The average baryon accretion rate is given by 



ab{t)^i}b 



ml 



SttG, 

= 1.2/i3Af0yr~^Mpc"^ 



dfs\ 



dz 

J}b_ 
0.044 



x(l + 2)v/r2A + 0,„(l + z)3 



dfs^ 



dz 



(26) 



where /struct is the fraction of total mass in star forming 
halos. It is given by 



J^^^^^JMMfps{M,z) 
/o"dAfA'//ps(A/,z) ' 



(27) 



and /ps is the Press-Schechter mass function. 

Outflows from the ISM into the IGM are described 
using the following relation 



oit) 



dm(j>{m) ■ ^[t{z) 



2e , 

i)\ ■ Ey,n-,(m), 
(28) 

where, as before, 4> is the stellar IMF, and i?kin is the 
kinetic energy released by the explosion of a star of mass 
m. The quantity Dose is the "average escape velocity," 
which is calculated according to 



J^^^^ dMfpsjM, z)M{2GM/R, 
Jm . dMfpsiM,z)M 



(29) 



The ejection of enriched gas from stars into the ISM is 
given by 



eit) 



dm (/)(to) ■^'[t(z) — r(m)] • [m~mr(m)], (30) 



as above, where TOu is the upper limit on mass of stars 
that explode and produce supernova remnants, and 
is the remnant mass. 

2.5. Ionization and thermal evolution of the 
inhomogeneous IGM 

Our model for reionization and thermal history of 
the average IGM is essentially that developed by 
Choudhury & Ferrara (2005), with the main difference 
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being in implementation of sources. In this model, which 
matches a wide range of observational constraints, the 
IGM is gradually reionized and reheated by star-forming 
galaxies between z ^ 20 and z 6. We summarise the 
main elements of the model here and refer the reader to 
that paper for details. 

The model accounts for IGM inhomogeneities by 
adopting a lognormal density distribution with the evo- 
lution of volume filling factor of ionised hydrogen (H ii) 
regions Qmiiz) being calculated according to the method 
outlined by Miralda-Escude et al. (2000); reionization is 
said to be complete once all low-density regions (say, 
with overdensities A < Acrit ^ 60) are ionised. We 
refer to the redshift at which this happens as the red- 
shift of reionization, Zjeion- We follow the ionization and 
thermal histories of neutral and H ii regions simultane- 
ously and self-consistently, treating the IGM as a multi- 
phase medium. We assume that helium is singly-ionized 
together with hydrogen and contributes to the photoion- 
ization heating. We do not consider the double ionization 
of helium, which is thought to have occurred at z ~ 3 
(e.g., Worseck et al. 2011). 

The rate of ionising photons injected by the galaxies 
into the IGM per unit time per unit volume at redshift z 
is denoted by fipY^lz), which is determined by the star for- 
mation rate. Using our star formation model described 
above, we can write the rate of emission of ionising pho- 
tons per unit time per unit volume per unit frequency 
range, n^{z), as 

M^) = /esc [<(z^)*"(z) + iV™(i.)vI'™(z)] (31) 

where /osc is the escape fraction of UV photons, and 
A^^^(i/) and N}^^^{i') are the total number of ionising pho- 
tons emitted per unit frequency range per unit stellar 
mass from Population I/II and Population III star form- 
ing halos respectively. These quantities are calculated by 
integrating the stellar spectra over an appropriate range. 
5'^^(z) and ^'^^'(z) are the cosmic Population I/II and 
Population III star formation rates. At any redshift. 



*(z) = «'"(z) + ^-"1(2). 



(32) 



The total rate of emission of ionising photons per unit 
time per unit volume, nph(z), is obtained by simply inte- 
grating Equation (31) over the correct frequency range. 
The hydrogen photoionization rate is then given by 



TMz) = {l + zf 



dv\{z,v)'hi,{z)a{v)^ (33) 



where is the hydrogen photoionization cross-section, 
and A(z, v) is the redshift-dependent mean free path of 
UV photons. The mean free path is governed by the 
distribution of high density neutral clumps and repro- 
duces the observed density of Lyman-limit systems at 
low redshift. The photoionization rate is used with the 
average radiative transfer equation to calculate the evo- 
lution of Qmi{z) in the pre-reionization universe, and 
that of Acrit (-z) in the post-reionization universe. Note 
that this approach assumes ionization equilibrium and 
local absorption of UV photons. 

The mean free path is defined as the average distance 
between high density regions of the inhomogeneous IGM. 



It can be written as (Choudhury & Ferrara 2005) 

Ao 



A(z) = 



[l-Fy(Aerit)]2/3' 



(34) 



where (Acrit) is the volume fraction of the IGM occu- 
pied by regions with density less than Acrit, given by 



-Fy(Ac 



dAP(A), 



(35) 



where -P(A) is the probability distribution function 
of IGM overdensities, taken to have a lognormal 
form (Miralda-Escude et al. 2000; Choudhury & Ferrara 
2005). We take Aq to be proportional to the Jeans length 
(Equation 40), 

Ao = KiipflTj, (36) 

and constrain the constant of proportionality, Anifp,o- us- 
ing the number of Lyman-limit systems at low redshift, 
which is related to the mean free path by 



^{j:)X{z)H{z){l + z) 



(37) 



We assume that the mean free path is independent of 
the photon energy up to the ionization threshold of He li 
and is proportional to (v /vnciiY'^ at higher energies 
(Choudhury & Ferrara 2005). 
The IGM temperature evolves as 



dT 
'dt 



-2H{z)T 



T dX 



HI 



1 + Xm dt 



dE 



3fcBnf,(l -|- z)^ dt ' 



(38) 



where dE/dt is the total photoheating rate and Ahi is 
the H I fraction. The photoheating term accounts for 
heating by the UV background and cooling due to re- 
combinations and Compton scattering of CMB photons. 
The IGM temperature governs the thermal feedback on 
galaxy formation by stopping gas inflow on halos with 
mass less than Mmin(z). We take this to be the Jeans 
mass, given by 

4 



where rj is the comoving Jeans length 



2kBTj 



1/2 



(39) 



(40) 



where 7 is the adiabatic index and mp is the mass of 
the proton. Note that before reionization is complete, 
the H II regions are distinct. During this epoch, the 
temperature in H II regions is different from that in H i 
regions. The average temperature is then given by the 
mean of these two temperatures, weighted by Qmiiz). 
Similarly, both regions have different Mj^in- This is rel- 
evant in understanding the role that Afmin plays in our 
star formation prescription above. 

Note that in all of above, we evolve two distinct sets of 
haloes corresponding to haloes forming in H i and H 11 
regions respectively. Each set of haloes is affected by 
photoionization feedback differently as the temperatures 
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Table 2 ^^-i 
Free parameters in our model 

^ ry-2 



Parameter 


Fiducial value 


Constraint 


/esc 


0.04 


Te and Phi 


/* (Pop. II) 


0.02 


Te and Phi 


/. (Pop. Ill) 


0.04 


Te and Phi 


•^mfp,0 


1.7 X 10-3 




e 


1.0 X lO-'^ 


/b, struct (^ = 0) 


/gas. in 


0.7 


Bouche et al. (2010) 


C 


Eq. (17) 


Krumholz & Dekel (2012) 




10-4 Zq 





of H I and H ii regions are different. Any average prop- 
erty, such as the cosmic SFR, at a given redshift is an 
average of these two sets weighted by the volume fill- 
ing fraction of H i and H ii regions. This distinction is 
unimportant at redshift z < Zrcion, when there are no H i 
regions anymore. 

Given the above model, we obtain best-fit parame- 
ters by comparing with various observations such as 
the redshift evolution of photoionization rate obtained 
from the Lya forest (Bolton & Haehnelt 2007), the 
electron scattering optical depth (Larson et al. 2010), 
and the total baryon fraction in structures at z = 
(Fukugita & Peebles 2004). These parameters are sum- 
marised in Table 2 along with their fiducial values and 
the constrains used to obtain them. 

3. RESULTS 

3.1. Global properties of the chemical evolution models 

Our model for galaxy and IGM evolution is constrained 
to produce the observed star formation rate (SFR) den- 
sity evolution, the fraction of total baryon density in 
haloes at low redshift, and the reionization history of 
the IGM as measured by the Thomson scattering optical 
depth to the last scattering surface (re) and the hydrogen 
photoionization rate (Fhi). Figure 1 shows the evolution 
of the cosmic SFR density. The solid black curve shows 
the prediction of our model and the magenta data points 
show observational measurements from a compilation by 
Hopkins & Beacom (2006). All three models reproduce 
the observed rise in the SFR between z = and z ^ 3. 
The total SFR prediction of the three models are virtu- 
ally identical. As we discuss in detail below, this is be- 
cause the contribution from Population III stars is very 
small. However, given the error bars on the measure- 
ments and the simplicity of our model, we consider the 
match to be adequate. Note that recent measurements 
of the SFR density at z - 7-10 by Bouwens et al. (2010) 
and Gonzalez et al. (2011) are smaller than our model 
prediction at these redshifts. However, this could pos- 
sibly be due to the detection limit in the observations, 
which miss the large contribution from faint galaxies. 

The star formation history shown in Figure 1 is also 
consistent with the measurement at z = of the stellar- 
to-halo mass (SHM) relation by Moster et al. (2012). 
The model prediction is in reasonable agreement with 
the the low-mass slope, characteristic mass (M = 10^^ 
M0), and normalisation of the measured stellar-to-halo- 
mass relation. The match at the high mass end is some- 
what poor, primarily due to lack of AGN feedback in 
our model. In massive haloes, AGN feedback can heat 
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Figure 1. Evolution of the eosmic star formation rate (SFR) den- 
sity in our models. Magenta points are observational measurements 
from a eompilation by Hopkins & Beacom (2006). In the bottom 
panel, the solid black curve shows the total (Pop. II -|- Pop. Ill) 
SFR density, defined in Equation (32). The two dashed curves 
show contribution of Pop. Ill SFR in two of our models: black 
and red curves show models 1 (1-100 Mq), and 3 (100-300 Mq) 
respectively. The top panel shows the fractional contribution of 
Pop. Ill SFR to the total SFR. 

the gas in the halo and prevent cooling and star forma- 
tion. We ignore this inconsistency since we are mainly 
interested in low mass galaxies in this work. Our fit to 
the SHM relation also moderately worsens at high red- 
shift. This agrees with the results of other semi-analytic 
models of galaxy formation, which also find it difficult to 
quantitatively reproduce the SHM relation (Sawala et al. 
2011; Guo et al. 2010) at both the high mass as well as 
low mass ends. However, there are several uncertainties 
in the high-redshift measurement of the SHM relation it- 
self, including those in the stellar mass estimation, faint- 
end luminosity function, and the halo mass function. As 
a result, we treat the SHM relation as only a qualitative 
check of the results of our models. 

Figure 2 shows the thermal evolution of the IGM in our 
models. The solid curve in left panel shows the evolu- 
tion of the temperature of the IGM at mean density (Tg). 
This is the temperature that is usually measured from ab- 
sorption spectra of high redshift quasars (Schaye et al. 
2000; Caluraetal. 2012). At redshift z > 20, Tq de- 
creases adiabatically oc (1 + z)^. However, as the HII 
regions around first galaxies form, this temperature grad- 
ually rises to about 2 x 10^ K at z = Zroion ~ 6. Once 
reionization is complete, only the highest density regions 
in the IGM are left to be ionized. Thus the mean den- 
sity regions evolve unaffected and Tq then decreases to 
~ 10^ K. At redshift higher than Zreiom is the mean 
of temperatures in H i and H ii regions, weighted by the 
volume filling factor. As a result, at these redshifts tem- 
peratures in these regions arc more relevant. In Figure 
2 (left panel), the red dashed curve shows the tempera- 
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Figure 2. Evolution of IGM temperature (left panel) and the filterini 
panels, the blue dashed curve shows the evolution in H I regions and the 
is shown by the solid black curve in both panels. 

ture evolution in H ii regions, and the blue dashed curve 
shows the same in H i regions. The average H ii region 
temperature climbs to a high value of a few time 10^ K 
at z ~ 20 when the first sources are turned on. How- 
ever, this does not affect Tq at these redshifts as the H ii 
region volume filling factor is small. The average tem- 
perature in H I regions continues to decrease with time 
due to adiabatic cooling. 

The thermal evolution of the IGM dictates the evolu- 
tion of Minin, the minimum mass of star- forming haloes, 
which is shown in the right panel in Figure 2. This quan- 
tity implements photoionization feedback in our model. 
As described in the previous section, we calculate Afmin 
using a Jeans criterion. At a given redshift the minimum 
mass of star-forming haloes is either given by the atomic 
cooling threshold of 10''^ K, or by the local Jeans mass, 
whichever is higher. Other models of Mmin are present 
in the literature, which take into account the full thermal 
history of the IGM instead of the instantaneous temper- 
ature (Gnedin 2000; Okamoto et al. 2008). We ignore 
these improvements and focus on the Jeans method for 
simplicity; this approximation does not afi^ect our results 
strongly. The solid curve in Figure 2 (left panel) shows 
the average M^i^ evolution. It increases monotonically 
with decreasing redshift and is ~ 2x10^ M© at z 



mass (right panel) in our models, shown here for model 3. In both 
red dashed curve shows the evolution in H II regions. The average 



It is ~ 6 X 10* Mq at redshift 1. It is important to note 
that this is only the average Mmin; the minimum mass is 
different in H i and H li regions. These are shown by the 
blue and red curves in Figure 2. Thus, galaxies forming 
in these regions are affected differently. The minimum 
mass in H ii region grows to about 10^ M© as soon as 
the first sources are turned on and the temperature of 
these regions is boosted to ^ 10* K. On the other hand, 
the minimum mass in H i regions increases much more 
slowly with decreasing redshift as it is always determined 
by the atomic cooling threshold alone. The relative con- 
tribution of the two regions to the solid black curve in 




Figure 3. Hydrogen photoionization rate in our models (solid 
curve), compared with measurements from the Lya forest. Blue 
data points show measurements by Meiksin & White (2004), red 
points arc by Bolton &; Haehnelt (2007), and magenta points are 
by Faucher-Giguere et al. (2008). 



Figure 2 is determined by the filling factor of H ii re- 
gions. Of course, this distinction is not important in the 
post-reionization phase, wheren there are no H i regions. 

As seen in the thermal evolution above in Figure 2, 
IGM reionization in our model is gradual. In this pic- 
ture, reionization begins at z ~ 20 and is 90% complete 
by z 7. This gradual change in the ionization state 
of the IGM helps us simultaneously reproduce the ob- 
served Thomson scattering optical depth to the last scat- 
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Figure 4. The mass-mctallicity relation prediction from our mod- 
els (solid curves) compared with observational measurements (data 
points). Measurements by Erb et al. (2006) at 2 = 2.27 are shown 
in black and those by Maiolino et al. (2008) at z = 3.7 arc shown in 
red. Note that the agreement between our model and the observa- 
tions here is only between their slopes, not between the amplitude, 
which are uncertain and have been rescaled as described in the 
text. 

tering surface Te = 0.089 and the measurements of the 
hydrogen photoionization rate from observations of the 
Lya forest (Meiksin & White 2004; Bolton & Haehnelt 
2007; Faucher-Giguere et al. 2008). The hydrogen pho- 
toionization rate in our model 3 is shown in Figure 3; 
our other models (1 and 2) give very similar results. 
The model prediction match reasonably well with the 
observational measurements. The photoionization rate 
increases rapidly as UV photon sources build up at z > 
•Zroion- Thcrc is a sudden jump at z = Zrcion when differ- 
ent H II regions overlap. This is because at this redshift, 
a given point in the IGM starts "seeing" multiple sources, 
which rapidly enhances the UV photon mean free path, 
thereby affecting the photoionization rate. Thus, our 
IGM reionization model simultaneously reproduces the 
two most robust observational constraints on reioniza- 
tion: (1) the Thomson scattering optical depth r to the 
last scattering surface, and (2) the hydrogen photoion- 
izaton rate, Fhi. Addition of a second stellar population 
at high redshift allowed us to produce the correct value 
of T without increasing Fhi at low redshift. Also, the 
underlying star formation physics in our model causes 
the time scale of star formation in a halo of given mass 
to evolve differently from previous models. This has a 
significant effect on the evolution of the photoionization 
rate (Chin & Ostriker 2000; Samui et al. 2007). 

We now turn to the mass-metallicity relation predic- 
tion in our models, which is shown in Figure 4. We find 
that the model predictions match quite well with ob- 
servations at z = 2.27 (Erb et al. 2006) and z = 3.7 
(Maiolino et al. 2008). The oxygen abundance increases 
with stellar mass for galaxies in the intermediate mass 
range. At the high mass end (M* > IO^^Mq) we find 
a slow decline in the oxygen abundance. This is be- 



Table 3 

Logarithmic yields of different species in 
our three models, obtained by using 
Equation (41). 



Species Model 1 Model 2 Model 3 



c 


-2.34 


-1.72 


-1.51 


N 


-5.42 


-6.12 


-6.49 


O 


-1.59 


-0.94 


-0.58 


Si 


-3.04 


-3.17 


-1.29 


Fe 


-2.71 


-4.31 


-1.87 


Zn 


-5.95 


-7.16 


-6.98 



cause the gas accretion rate in these haloes is high, 
which results in dilution of their ISMs. Also, it is dif- 
ficult for outflows to remove gas from the deeper poten- 
tial wells of these halos. Note that the measurements 
of the mass-metallicity relationship are quite uncer- 
tain (Kewley & Elhson 2008; Krumholz & Dekel 2012). 
Therefore, to compare our estimate of metallicity to the 
observations we force the zero point of the predictions to 
match with that of the observations at z = 2.27. The 
model mass-metallicity relationship has a slightly differ- 
ent slope than the observational points at z = 2.2. How- 
ever, given the simplicity of our model, and our focus on 
ratios of abundances in this work, we take the level of 
agreement to be sufficient. We will also see below that 
when the average ISM metallicity of galaxies is averaged 
agrees perfectly with the mean ISM metallicity inferred 
from observations of damped Lya absorbers (DLAs) . 

In summary, our average galaxy formation and chemi- 
cal evolution model is consistent with observable quanti- 
ties, including the observed Thomson scattering optical 
depth to the last scattering surface, the observed hydro- 
gen photoionization rate inferred from the Lya forest, the 
slope of the mass-metallicity relation of galaxies at vari- 
ous redshifts, and the evolution of cosmic star formation 
rate density. It also agrees well with the low-mass slope, 
characteristic mass, and normalisation of the stellar-to- 
halo-mass relation, and the evolution of the mean ISM 
metallicity inferred from observations of DLAs. 

3.2. Effect of the Population III stellar IMF 

To study how the Population III stellar IMF affects 
above results, we now consider three different Popula- 
tion III stellar IMFs shown in Table 1. In each case, 
our Population II IMF is fixed to 0.1-100 Mq Salpeter. 
Each of these three models is calibrated independently, 
as described above. Our goal is to understand how the 
difference in these Population III scenarios is reflected 
in chemical evolution of galaxies. The 1-100 Mq and 
100-300 Mq IMFs are selected to represent two extreme 
possibilities: if metal-free gas can fragment due to dy- 
namical effects the IMF will favour small stellar masses, 
but in the absence of fragmentation, Population III stars 
could have masses of the order of ~ 100 Mq. The inter- 
mediate (35-100 M0) IMF is chosen to consider the effect 
of core collapse supernovae. If the distinct chemical sig- 
natures of these three IMFs turn out to be observable in 
high redshift galaxies, then such observations could act 
as a probe of the Population III IMF. 

The contribution of Population III stars to the total 
cosmic SFR density is shown in Figure 1, in which the 
bottom panel shows the evolution of the Population III 
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halo mass [N/0] 

Figure 5. Abundance ratios and DLA distribution functions predicted by our model at ^ = 6. Black, blue and red curves correspond to 
models 1,2, and 3 respectively. 




Figure 6. Abundance ratios and DLA distribution functions predicted by our model at 2 = 2, shown on the same scale as Figure 5. 
Black, blue and red curves correspond to models 1, 2, and 3 respectively. 
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SFR (\E'm(z)) and the top panel shown the fractional 
contribution to the total SFR (*(z)/*"i(z)). Evidently, 
Population III stars contribute very little in terms of 
SFR. Furthermore, any Population III contribution is 
essentially zero by redshift z ~ 7. The initial burst 
of Population III stars is sufficient to enhance the ISM 
metallicity beyond Zf-rit in any halo mass bin. When this 
happens, the corresponding galaxy stops forming Popu- 
lation III stars. Since most of the mass in the universe 
is contained in Af» haloes, the globally averaged Popula- 
tion III SFR starts declining as soon as these haloes cross 
the Zcrit threshold. Since the IMF in model 1 has a lower 
metal yield than that in model 3, haloes take longer to 
cross the Zait threshold. Therefore, the Population III 
SFR contribution is larger in model 1. These results are 
in good agreement with the hydrodynamic simulations 
of Wise et al. (2012), who use only very massive Popu- 
lation III stars. We find that while the amplitude of the 
Population III SFR is sensitive to this effect, the lowest 
redshift with nonzero cosmic Population III SFR den- 
sity is more dependent on the evolution of the minimum 
mass of star-forming haloes, M^un- Further, the small 
contribution of Population III stars to the global SFR 
is reflected in the fact that the effect of changing the 
Population III stellar IMF on the reionization and ther- 
mal history of the IGM is negligible in our model. As 
a result, the temperature evolution, as well as the pho- 
toionization rate evolution is practically independent of 
the Population III IMF. 

However, we find that although the efi'ect of Popula- 
tion III on the SFR and the IGM reionization and ther- 
mal history is small, it has a significant efi'ect on the 
chemical properties of low mass galaxies. This is seen 
in Figure 5 (left column), which shows abundance ratios 
of various metal species as a function of halo mass at 
z = 6 for our three models. We focus on [C/Fe], [0/Si], 
[Zn/Fe], and [N/0]. It is seen that in haloes with mass 
greater than a few times 10^ Mq, all four abundance 
ratios are constant, and equal to the Population II val- 
ues. Changing the Population III IMF has no effect on 
the abundance ratios in these haloes. This is because in 
these haloes. Population III star formation has ceased at 
a much earlier time and the subsequent Population II star 
formation has wiped out any chemical signature of Pop- 
ulation III star formation. The constancy of the abun- 
dance ratios is due to the fact that for a time-independent 
IMF, abundance ratios are equal to the ratio of corre- 
sponding chemical yields (Tinsley 1980; Pagel 2009). In 
low mass haloes, on the other hand, chemical signatures 
of Population III star formation have not yet been wiped 
out. As a result, abundance ratios in these haloes are 
different in each of our Population III models. Thus, the 
chemical signature of Population III stars slowly emerges 
as we move towards low halo masses. Our model suggests 
that this is probably the best hope of constraining the 
Population III SFR density on cosmic time scale and the 
Population III IMF. 

When we change the Population III model, the abun- 
dance ratios in low mass haloes change significantly. For 
example, the [0/Si] ratio shows a change of more than 
1 dex when we change from model 2 to model 3. This 
can be seen in Table 3, which shows the Population III 
yields of O and Si in all three models. Here, the yield of 
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Figure 7. Comparison of ciiemical evolution in haloes with three 
different masses, shown here for model 3. The dashed line in the top 
panel shows the If II region filtering mass from Figure 2, whereas 
the dotted line shows the evolution of the average filtering mass 
seale. The black, red, and blue eurves denote three different haloes. 
The dashed vertical line marks z = 6, which corresponds to the 
redshift in Figure 5. The spike in the plot is a result of the most 
massive stars in the model. 

species i is defined as 

J dm ■ <j){m) ■ m ■ p,{m) 

= — 'n 7r~\ — ' (^^^ 

J am ■ m ■ (jjim) 

where pi (m) is the fraction of initial stellar mass m that is 
converted to species i. (Table 3 shows log(yi).) The ratio 
of oxygen yield to silicon yield is comparable to solar in 
model 1. However, it is much lower than solar in model 3, 
and much higher than solar in model 2. This is because 
high mass Population HI in model 3 stars produce large 
amounts of silicon during the 0-burning process. This 
is reflected in the [0/Si] ratio shown in Figure 5. Sim- 
ilar considerations explain the differences seen in other 
abundance ratios in low mass haloes in the three models. 
As we move towards lower redshift, we would expect to 
see the abundance ratio values to move away from their 
Population HI values towards the Population II values 
even in small mass haloes. This is seen in Figure 6 (left 
column) which shows the same abundance ratios as Fig- 
ure 5 at z = 2. The abundance ratio values in low mass 
haloes are much closer to the Population II values. 

3.3. Why do low mass galaxies retain an imprint of 
Population III stars? 

As described above, our model predicts that Popula- 
tion HI stars will leave an imprint on the ratio of abun- 
dances of metals in low mass haloes, while any such im- 
print is wiped out in high mass haloes by subsequent 
Population II star formation. In this section, we dis- 
cuss the reason behind this contrast between low mass 
and high mass haloes. In a nutshell, it is caused by a 
difference in Population Ill-to-II transition redshifts of 
these haloes, which in turn is a result of photoionization 
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feedback. 

The solid curves in Figure 7 (top panel) show the aver- 
age growth of three different haloes from H ii regions in 
one of our models (model 3). These haloes have masses 
of lO^-^ Mq, 10^°-5 Mq, and 10"-^ Mq respectively at 
z = 1. Their growth is described by Equation 1. The 
dashed line in this panel shows the evolution of the fil- 
tering mass Afmin in H il regions in this model. It can 
be seen that the three haloes satisfy the star formation 
criterion, M > Afmin at three different redshifts, approx- 
imately z = 10, 20, and 25 respectively. Thus each halo 
forms its first stars, which are Population III, at very dif- 
ferent times. This has an effect on the chemical evolution 
of these three mass bins, which is shown in the bottom 
panel of Figure 7. In this panel, the three solid curves 
show the evolution of [0/Si] in the same halo mass bins 
as the top panel. (Wc still focus on model 3.) Since the 
Population III IMF in model 3 (100-300 Mq) produces 
negative values of [0/Si], we expect a negative value of 
[0/Si] in each halo mass bin when Population III stars 
form. However, as time progresses. Population II star 
formation takes over in each mass bin, because of which 
the value of [0/Si] approaches -1-0.2. This is seen to hap- 
pen in all three halo mass bins. However, this happens at 
different redshifts for different halo masses because they 
have different redshifts of Population HI star formation 
(and therefore also the redshifts of Population Ill-to-II 
transition) . The redshift at which Population HI star for- 
mation first occurs is the redshift at which the condition 
M > Mnun is first met, and is thus governed by the ther- 
mal feedback. As a result, at a given redshift, low mass 
haloes retain memory of their Population HI enrichment 
history. This is highlighted by the vertical line marking 
z = 6 on the bottom panel of Figure 7. It can be seen 
that at this redshift, the two higher halo mass bins have 
already converged to [0/Si] = 0.2, but the lower mass bin 
still has [0/Si]= —0.3. It is this effect that causes the 
trends seen in Figure 6. 

We now sec why changing the Population HI stellar 
IMF will selectively affect low mass haloes. Figure 8 
shows the evolution of [0/Si] for the same three haloes 
as Figure 7 in each of our three Population HI models. 
As expected from the yields in Table 3, [0/Si] for Pop- 
ulation HI stars in models 1 and 2 is greater than 0.2, 
where that for model 3 is negative. The curves for mod- 
els 1 and 2 show an enhancement in [0/Si] during the 
first star formation episode in each halo mass bin. But 
this occurs at different redshifts for different halo masses, 
and explains the variation in the trends seen in Figure 
6. An alternate way of understanding this is that in a 
closed system with a non-evolving IMF, all abundance 
ratios get locked at the value given by the ratio of their 
respective yields (Tinsley 1980; Pagel 2009) 



where the yields arc defined in Equation (41). While 
this is strictly valid only under the instantaneous recy- 
cling approximation, it holds for real galaxies to a very 
good approximation (Tinsley 1980). As a result, abun- 
dance ratios will change only when there is a change in 
the IMF. This is exactly what happens when the stellar 
IMF in a galaxy changes from Population HI to Popula- 
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Figure 8. Evolution of [O/Si] in our three models in tiie same 
three mass bins as in Figure 7, shown here in the same colours. 
The dotted, dot-dashed, and solid curves show evolution in models 
1,2, and 3 respectively. 

tion II. For low galaxies, this happens at a lower redshift. 
Therefore, they retain the Population HI signatures at 
these redshifts. 

It is often mentioned in the literature that an "odd- 
even" pattern in the abundances is a signature of Popu- 
lation HI stars. This refers to the fact that in metal-free 
stars the production of species with odd nuclear charge 
is preferentially suppressed relative to their solar values. 
As we have seen above, the low mass haloes in our mod- 
els retain memory of Population HI star formation and 
therefore also show the odd-even pattern in elemental 
abundances. However, it is important to understand that 
odd-even pattern is not unique to Population HI and also 
occurs in metal-poor Population II stars (Cooke et al. 
2011). Hence the crucial point here is that one really 
wants to search for relative abundance trends which differ 
from the Population II predictions, and this information 
is encoded in low-mass halos. 

3.4. Role of critical metallicity 

In all of above, the critical metallicity ^crit, at which 
a Population HI IMF changes over to a Population II 
IMF, plays a crucial role. However, the value of 
Zciit is currently a topic of debate (Klessen et al. 2012; 
Dopcke et al. 2012). When metallicity Z is greater 
than Zcrit, the availability of metal-line cooling and/or 
dust cooling is expected to lead to fragmentation in 
the ISM that could help form small mass stars with 
a Population II type IMF. In this paper, we have as- 
sumed Zcrit = W~^Zq, a value motivated by most stud- 
ies of fragmentation in metal-poor gas (Bromm et al. 
2001; Bromm & Loeb 2003; Frebel et al. 2007). How- 
ever, some studies have argued for a much smaller value 
of Zcrit = 10~^Zq (Klessen et al. 2012). Since this num- 
ber is not very well constrained, we comment on the effect 
on our results of variation in its value. 

We find that values of ZcHt smaller than IO^'^Zq have a 
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Figure 9. Evolution of dN/dX in our model. The red data points 
from measurements by Prochaska et al. (2005). 

very minor influence on our results. This is because in all 
three of our models, the metallicity of ISM in all haloes 
crosses Zcrit in a very short time after the first burst 
of Population III star formation. We explicitly checked 
this in our calculation by changing Zcrit to IO^^Zq and 
IQ^^Zq. In our models 2 and 3, in which the Popula- 
tion III metal yield is higher, there is no effect of a change 
in the critical metallicity. In model 1, lowering the criti- 
cal metallicity to 1Q~^Zq has an effect of decreasing the 
average Population III SFR density, but only for z > 20. 
When Zcrit is lowered, haloes in this model take shorter 
time to stop forming Population III stars. These results 
arc in qualitative agreement with those of Salvadori et al. 
(2007) , who also found that the influence of Zcrit depends 
on the Population III stellar IMF. Given the negligible 
effect of changing Zcrit, we keep Zcrit = IO^^Zq in this 
paper. 

3.5. Observability as DLAs 

As we saw above, low mass galaxies (halo masses ^ 10^ 
Mq) at high redshift (z ~ 6) are likely to be useful in con- 
straining Population III IMF and SFR. However, these 
galaxies are difficult to observe at high redshift. Even if 
we could observe them, it would be very difficult to mea- 
sure their metallicity via line emissions thanks to highly 
uncertain modelling of H ii regions. As a result, we need 
to look for an alternative way of observing these galax- 
ies and measuring their metallicities. We now consider 
their observability as absorption systems in high-redshift 
quasar spectra. 

DLAs comprise of neutral hydrogen reservoirs at high 
redshift. Measurements of their metal content is useful 
in understanding chemical evolution of galaxies. Metal 
abundances of DLAs have been measured up to z = 6. 
Using this data, an "age-metallicity relation" can be con- 
structed for galaxies. Moreover, chemical abundance 
measurements in DLAs are very accurate. This is be- 
cause in the cold H i reservoirs of DLAs, most metal 
species are in their singly ionised state, which can be 
probed in the optical from ground-based telescopes. Er- 
rors in DLA H i column density measurements are typi- 
cally quite low, around 0.05 dex. Therefore errors in the 
corresponding metal abundance measurements ([M/H]) 
is also low, around 0.1 dex. Thus, observations of 



DLAs can be used to measure gas-phase metallicities 
at large cosmological lookback times with high precision 
(Wolfe et al. 2005). Contrastingly, abundance measure- 
ments in galaxies in the high redshift universe is quite 
uncertain because modeling of ionized gas giving rise to 
line emission in H ii regions is complicated and involves 
many uncertainties, particularly at high-redshift where 
one only covers a handful of rest-frame optical lines in the 
observed- frame near-IR (Kewley & Ellison 2008). Fur- 
thermore, in DLAs, relative abundances can still be mea- 
sured accurately deep into the reionization epoch (z > 6) 
using transitions red ward of Lya, even though Gunn- 
Peterson absorption precludes measurement of neutral 
hydrogen (Becker et al. 2012). 

For any measurable property p of DLAs at high red- 
shift, we can calculate the line density distribution of 
DLAs as a function of that property. This tells us the 
contribution of systems with different values of the prop- 
erty p to the total sample of DLAs. If we assume a 
one-to-one correspondence between a halo's total mass 
and its H i absorption cross section ("size"), the DLA 
line density distribution can be written as 

(PN , , , , dl dM , 
^.iV(A./).E(M).--.(l + z)3. (43) 

Here, p is the property in consideration (e.g., metallicity), 
dl = cdt is a length element, and dX is an absorption 
length element given by 



dl 



dX Hail + zf 



(44) 



The halo mass is given by M, and N{M) is the comoving 
number density of halos, i.e., the halo mass function. 
Wc then assign an absorption cross section to each halo, 
using the fitting function 



I](M) = So 



M 



1 



M 
Mo 



a-2 



(45) 



where the constants take the values of a = 0.2, Mq ~ 10^ 
Mq, and So = 40 kpc at z = 3 (Pontzen et al. 2008; 
Font-Ribera et al. 2012). This function provides a good 
fit to the absorption systems in hydrodynamical simula- 
tions of Pontzen et al. (2008) at z = 3. (We had to as- 
sume a smaller value of a to get a good match to observed 
dN/dX.) At other redshifts, we calculate the values of 
these constants by assuming that haloes with same cir- 
cular velocity at any redshift have the same absorption 
cross section (Font-Ribera et al. 2012). Further, since 
properties like gas mass, stellar mass, and metallicity 
also have one-to-one correspondence with halo mass in 
our model, we can easily calculate dM/dp. 

Once the absorption cross-section is defined, we can 
look at the average evolution of ISM metallicities. The 
average ISM metallicity at a given redshift is 



Z(z) 



/ dm ■ N{m, z) ■ ZisM{m, z) ■ Y.{m, z) 
/ dm ■ N{m, z) ■ S(m, z) 



(46) 



Figure 10 shows this "age-metallicity" relation in our 
model in comparison with metallicity estimates for 
100 damped Lya systems (DLAs) with z < 4 by 
Prochaska et al. (2003). Our models arc in good agree- 
ment with the mean metallicities of the observed DLAs. 
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Figure 10. The "age-metallicity" relation in our model, which 
is given by an absorption cross section weighted average of gas 
phase metallicity of all mass bins. Red crosses are measurements 
by Prochaska et al. (2003). 

The integral of Equation (43) over all values of p is 
just the total line density of DLAs, dN/dX. The evolu- 
tion of dN/dX in our model is shown in Figure 9. Our 
results are in reasonable agreement with the values ob- 
tained from the Sloan Digital Sky Survey (SDSS) DR5 
by Prochaska et al. (2005). 

We plot the line density distribution of DLAs as a func- 
tion of various metal abundance ratios in Figure 5 (right 
column). We focus on [C/Fe], [0/Si], [Zn/Fe], and [N/0] 
as before. As discussed before, all high mass haloes in 
our model have the same value for the abundance ratios, 
given by the ratio of yields for Population IL Therefore, 
when our Population III IMF is very close to the Popu- 
lation II IMF, as in model 1 (1-100 Mq), the DLA line 
density distribution is a straight line centred at the con- 
stant value of the abundance ratio. In other models too, 
such as models 2 (35-100 Mq) and 3 (100-300 Mq), the 
high mass haloes have the same values of abundance ra- 
tios because their Population III signatures have already 
been washed out by subsequent Population II star for- 
mation. As a result, these systems also fall on the same 
line. However, the low mass haloes in these models have 
different abundance ratios. Therefore, the line density 
distribution at low halo masses deviates from the verti- 
cal line, over a range of values of the abundance ratio. 
For example, the [0/Si] abundance ratio for model 3 (red 
curve) stays on the vertical line for high mass haloes but 
deviates from it over [0/Si] = -0.5 to [0/Si] = 0.2. The 
line density distribution for model 3, seen as a function 
of [0/Si], has a very different behaviour. In this model, 
low mass haloes have a negative value of [0/Si]. Hence, 
the DLA line density distribution is spread out over neg- 
ative values of [0/Si]. Also, we would expect that with 
time, this distribution approaches the vertical line, which 
corresponds to the Population II value. This is seen in 
Figure 6 (right column), which shows the line density 
distributions at z = 2. It seen that the tails of the line 
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Figure 11. Top panel shows the predicted distribution of [O/Si] 
values for Population II (black) and Population III (red; model 
2) at 2 = 6. (This distribution is taken from Figure 5.) The 
bottom panel shows a histogram of a random sample of 10 [O/Si] 
measurements taken from the Population III distribution in the 
top panel. An additional Gaussian error of 0.1 dex is added to 
each sampled point. A KS test shows that this sample rejects 
the Population II distribution at approximately Aa {D = 0.9, p = 
0.00017). This suggests that it is possible to detect the effect of 
Population III star formation even with a relatively small sample 
of DLA abundance ratio measurements with an accuracy of about 
0.1 dex. 

density distribution are significantly shorter at z = 2 as 
compared to z = 6, which is because by this time the 
chemical signature of Population HI star formation is 
wiped out even in most lower mass halos. 

Thus, the distribution of DLAs in the abundance-ratio 
space at sufficiently high redshift is sensitive to the Pop- 
ulation III IMF. In the absence of any change in stellar 
yields over time, this distribution will be a "delta func- 
tion" in the abundance-ratio space (assuming corrections 
for the impact of dust on relative abundances). However, 
the addition of a new stellar population with different 
chemical yields, spreads DLAs out in the abundance- 
ratio space. The form of this spread can constrain the 
IMF of Population HI stars, while its evolution can con- 
strain the Population HI SFR history. 

Further, to check whether this effect is detectable and 
whether it can indeed be used to probe the Population HI 
IMF, in Figure 11. we show a simulated data set of 10 
[O/Si] measurements at z = 6 taken from the Popula- 
tion HI distribution corresponding to model 2 (shown in 
the top panel of Figure 11 and also in Figure 5.) This 
data set is produced by random sampling the predicted 
distribution of [O/Si] values in the model, and adding a 
Gaussian error of 0.1 dex to each sample point. An error 
of 0.1 dex corresponds to the typical accuracy with which 
abundance ratios in DLAs are measured (Becker et al. 
2012). Using a KS test, we find that this sample re- 
jects the Population II distribution at 4a {D = 0.9, 
p = 0.00017). A set of 100 samples (each of size 10 
and a Gaussian error of 0.1 dex) prefered the Popula- 
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tion III IMF at 3.8cr on average {(D) = 0.8, (p) = 0.007). 
This suggests that the effects of Population III IMF on 
DLA abundance ratio is significant enough to be de- 
tectable with just 10 measurements accurate to about 
0.1 dex^. A sample of 20 measurements, also with accu- 
rancy of 0.1 dex resuh in (D) = 0.85, (p) = 10"^ (about 
5a). Note that this test using a single abundance ra- 
tio ([0/Si]). In practice, using multiple ratios as shown 
in Figure 5 can further improve the significance of the 
constraints. Also note that the metallicity of these low 
mass halos (M lO^M©) with Population III signatures 
is high enough to produce detectable lines at z ^ 6. 
For example, we can estimate the O i column density 
by (Salvadori & Ferrara 2012) 



3 Mo/fJ-om-p 



2lT 



H. I ,,,, 



(47) 



where //q is the atomic weight of oxygen and r = 
arvir is the gas radius, where usually a ^ 0.8 
(Salvadori & Ferrara 2012). The oxygen mass Mq is cal- 
culated in each halo in our chemical evolution model. 
This gives Nqi > 10^^ cm~^ for a = 0.8. This value of 
the column density is already higher than that measured 
in several z = 6 DLAs by Becker et al. (2012), which 
span the range A^oi = lO^'^'^^-lO^^ '*'' cm~^. 

It is important to understand that the distribution of 
relative abundance values at high redshifts is used here 
and not the mean value at those redshifts. In fact, as 
we have seen, all distributions in Figure 5 arc anchored 
on the Population II values of the respective abundance 
ratios. This is because even at high redshift, high mass 
haloes have already moved to their Population II star 
forming stage; only low mass haloes are expected to carry 
the memory of their Population III enrichment. This can 
be understood by calculating the cross-section-weighted 
average value of different abundance ratios, given by 



(P) 



dXdp 



■ p ■ dp ■ dX, 



(48) 



where p =:[Mi/M2] is the abundance ratio of two species 
Ml and M2, and dX is related to the redshift bin dz over 
which the average is taken, given by Equation (44). Fig- 
ure 12 shows the average value of [C/Fe] and [0/Si] cal- 
culated in this fashion, corresponding to our three Pop- 
ulation III models. Observational measurements of sev- 
eral DLAs from Becker et al. (2012) are also shown. It is 
seen that the average values predicted by our model are 
consistent with the observations, regardless of the Pop- 
ulation III IMF used. (Note that the observational data 
shown here is uncorrected for effects of dust depletion.) 
The most striking feature in Figure 12 is the lack of evo- 
lution over the large redshift range from z > 6 to z ^ 2. 
This is explained in our model by the fact that all large 
haloes show constant abundance ratios, corresponding to 
Population II. A second feature of Figure 12 is that vary- 
ing the Population III IMF results in only a small change 
in the average value of the abundance ratio, as expected. 
It is only when we look at the distribution of these val- 

Any effects of dust depletion could in principle be higher than 
0.1 dex. However, the scatter in observations at low redshift {z ~ 2) 
is less than 0.2 dex (Becker et al. 2012), which suggests that our 
results are valid even when effects of dust are included. 
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Figure 12. Evolution of the average values of [C/Fc] and [O/Si] 
in our models. The black, blue, and red curves correspond to mod- 
els 1, 2, and 3 respectively. Data points in are from Becker et al. 
(2012) and references therein. It is seen that changing the Pop- 
ulation III IMF has little effect on the average values of relative 
abundances. However, the effect on the distribution of abundance 
ratio values is significant and detectable, as shown in Figure 11. 
This shows that it is the distribution of measurements that probes 
the stellar IMF and not the mean value of the measurements. 

ues at different redshift that the effect of Population III 
IMFs becomes visible. 

4. DISCUSSION 

We have shown that a consistent chemical evolution 
model that also incorporates global effects like reioniza- 
tion and photoionization feedback is consistent with a 
wide range of Population III scenarios. Furthermore, this 
approach uncovers patterns in the metallicity distribu- 
tion of high redshift galaxies which could potentially be 
useful in ruling out some Population III IMFs and SFRs. 

It is therefore important to understand the domain 
of validity of our model. Matteucci & Calura (2004) 
and Calura et al. (2009) have studied the evolution of 
the mass-metallicity relation and the chemical proper- 
ties of systems like dwarf galaxies using chemical evolu- 
tion modelling of the kind that we have used here. They 
find that if the instantaneous recycling approximation 
is avoided, an approach that correctly incorporates ef- 
fects of outfiows and inflows successfully reproduces all 
observed chemical properties of galaxies. De Lucia et al. 
(2004) and De Lucia & Borgani (2012) have also used a 
similar chemical evolution model in their semi-analytic 
model of galaxy formation that fits several key obser- 
vations very well. Furthermore, we focus on abundance 
ratios in this paper, which are more robustly predictable, 
given their simple dependance on the stellar IMF. Thus, 
given the uncertainties in high redshift observations and 
the simplicity of our model, the main conclusion of this 
paper regarding chemical signatures of Population III 
stars in high redshift galaxies is robust. 

A second uncertainty in our calculation is the pre- 
scription used for photoionization feedback. Firstly, 
we have used a simple Jeans prescription whereas it 
has been argued that more detailed prescriptions give a 
better fit to hydrodynamical simulations (Gnedin 2000; 
Okamoto et al. 2008). However, the difference between 
these feedback prescriptions is small, especially at z > 5 
(Barkana & Loeb 2001). A detailed understanding of ra- 
diative feedback in galaxy formation is yet to emerge. 
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primarily because of the large variety of approximations 
used in earlier work (Ciardi & Ferrara 2005). However 
the general consensus is that such a filtering scale in- 
creases from roughly lO*" Mq at z = 10 to 10^ Mq at red- 
shift z = 6. This is in good agreement with the values we 
obtain in our study. A tentative observational evidence 
of such effect is seen in "dark" dwarf galaxies in the vicin- 
ity of high redshift AGN (Cantalupo et al. 2012) and in 
the difference between the evolution of star formation 
rate in low mass galaxies and that in high mass galaxies 
(Weinmann et al. 2012). Some authors have used simple 
fitting functions for that agree quite well with our 

evolution of M^-m (Salvadori & Ferrara 2012). A similar 
approach has also been used in semi-analytical models of 
galaxy formation (e.g., Baugh 2006). It should also be 
noted that the evolution of M^i^ for z < 5 is not im- 
portant for our present study as by such low redshift the 
contribution of Population III stars to chemical evolution 
is negligible. 

Thirdly, our approach of evolving haloes and galaxies 
in the mean does not provide us with information re- 
garding the scatter around the distributions shown here. 
However, we note that semi-analytic models based on 
halo merger trees find quite small scatter around the 
mean distributions of mctallicities (around 1 dex at red- 
shift 0; De Lucia et al. 2004). Further, scatter does not 
affect abundance ratios; it can only affect the star forma- 
tion histories of individual haloes. Therefore the effect 
of scatter on our central prediction, that distribution of 
DLA abundance ratios probe Population HI IMFs, will 
be nearly zero. For example, positive values of [0/Si] at 
high redshift are inconsistent with Population HI stars 
in the PISN range. Observations of these systems would 
be a strong argument against high mass Population HI 
stars, regardless of the scatter in the observations, which 
will not affect the [0/Si] ratio. In fact, we want to high- 
light the fact that in this respect this approach of con- 
straining Population HI scenarios from observations of 
high redshift DLAs is possibly better than constraining 
them using observations of metal-poor stars in the Galac- 
tic halo. The mctallicity distribution of metal-poor stars 
is highly dependent on the peculiar assembly history of 
the Milky Way, which is not guaranteed to follow the "av- 
erage" assembly history that models require in order to 
constrain Population HI SFR and IMF. As such, metal- 
poor stars are useful only if they unambiguously show the 
signature of a pair instability supernova. However, even 
in this case they will only provide indication of a single 
Population HI star of a certain mass, not of the whole 
IMF or SFR.^ Similar arguments hold for observations 
of dwarf sphcroidals and ultra-faint dwarfs. 

We have assumed a constant UV photon escape frac- 
tion while calculating the IGM reionization history. This 
can potentially affect our conclusions if the escape frac- 
tion in Population HI star forming haloes is so large that 
reionization constraints begin restricting Population HI 
star formation parameter space. Indeed, it has been ar- 
gued that the escape fraction will have a strong increas- 
ing evolution with redshift (Kuhlen & Faucher-Giguere 
2012; Mitra et al. 2012). Also, the escape fraction could 

^ Another way in which mctal-poor stars can be of use is if we 
find a truly metal- free star. But this is difficult given the observa- 
tional limitations (Frebel & Norris 2011). 



potentially have different values for galaxies with Popula- 
tion II and Population HI star formation due to different 
spectral indices (Choudhury & Ferrara 2006). However, 
our assumption about the escape fraction is unlikely to 
affect our results significantly since we find that Popula- 
tion HI stars do not contribute much to reionization due 
to their early termination. Thus, the reionization history 
affects our conclusions only moderately. Note that the 
value of escape fraction used in our model is comparable 
to that deduced from luminosity function measurements 
in recent studies (Mitra et al. 2012). 

Finally, we comment about the role of dust. Dust af- 
fects the metal evolution of galaxies by preferentially 
depleting certain species like iron. We have not taken 
this effect into account. As such our results are only 
"production-side" estimates of the chemical abundance 
patterns. Dust obscuration is ignored in most DLA stud- 
ies (Wolfe et al. 2005). This difficult problem can be alle- 
viated by focusing on non-refractory species such as zinc 
and oxygen, which are not depleted on dust grains. We 
can also focus on "secondary elements" like nitrogen that 
depend on the square of the mctallicity, which increases 
their sensitivity to underlying chemical enrichment pat- 
terns. We defer a detailed study of the effect of dust 
depiction to a future work. 

5. CONCLUSIONS 
Our results are as follows: 

• We have developed a chemical evolution model of 
galaxies that is consistent with a variety of global 
observational constraints on galaxy and IGM evo- 
lution, such as the cosmic SFR density evolution, 
IGM thermal evolution, and hydrogen photoion- 
ization rate evolution. We calculate the minimum 
mass of star forming galaxies self-consistently. This 
model produces galaxies that lie on observational 
curves such as the stellar-to-halo mass relation at 
low redshift and the mass-metallicity relation. We 
then explore influence of Population HI stars on 
the predictions of this model, by varying the Pop- 
ulation HI stellar IMF. 

• We find that different Population HI stellar IMFs 
result in very different abundance ratio distribu- 
tions in low mass galaxies (halo mass < 10^ M©) 
at high redshift. This is because photoionization 
feedback suppresses star formation in these galax- 
ies till low enough redshift (z ^ 10), and memory of 
the initial generation of Population HI stars is re- 
tained. For some ratios, e.g., [0/Si], the variation 
is as large as 1 dex. This variation is much greater 
than the uncertainty in mctallicity measurements 
in DLAs. This effect is strong at redshift z >5 and 
grows weaker at lower redshift due to subsequent 
Population II star formation. 

• The influence of Population HI star formation is 
seen only in low-mass haloes. Changing the Popu- 
lation HI IMF has no effect on the abundance ratio 
in high mass haloes. This is as expected since in 
these haloes abundance ratios are determined by 
Population II stars alone. It is only as we move 
downwards in halo mass that the Population HI 
signature starts dominating. 



18 



Kulkarni ct al. 



• We modelled low-mass haloes as DLAs by assign- 
ing them a mass-dependent H i absorption cross- 
section. This model agrees with the observed 
mctallicity-rcdshift relation for DLAs and the ob- 
served DLA line density evolution as a function of 
redshift. Using this model, we predict the line den- 
sity distributions of DLAs as a function of different 
abundance ratios. We find that these distributions 
are anchored towards the Population II values of 
abundance ratios, but they show a significant de- 
pendence on the Population III IMF for z > 5. 
This dcpendance is reduced at lower redshift due 
to subsequent Population II star formation. 

• The dependance of DLA line density distribution 
on Population III suggests that the distribution of 
DLAs in abundance-ratio space at sufficiently high 
redshift can provide very good constrains on Pop- 
ulation III properties. The form of this distribu- 
tion can constrain the IMF of Population III stars, 
while its evolution can constrain the global Pop- 
ulation III SFR history. This constraint on the 
Population III could possibly be stronger than con- 
straints from other probes such as metal-poor stars 
and individual metal-poor DLAs. A simulated data 
set of just 10 DLAs at z ~ 6 measured with real- 
istic accuracy is able to constrain specific Popula- 
tion III IMFs with high confidence. This method 
of probing Population III stars does not rely on the 
measurement of the H i column density. It is there- 
fore useful at high redshift where Gunn-Pctcrson 
absorption precludes measurement of neutral hy- 
drogen. Relative abundances can still be measured 
accurately in this epoch using transitions rcdward 
of Lya. 

• Not only can the abundance ratio distributions dis- 
tinguish between high-mass and low-mass Popula- 
tion III IMFs, they can also discriminate between 
different low mass Population III IMFs. For ex- 
ample, the two low mass Pop III IMFs that we 
consider in this paper are clearly distinguished by 
their oxygen abundance patterns. 
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